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We evaluate the two-photon exchange (TPE) correction to the unpolarized elastic 
electron-proton scattering at small momentum transfer Q 2 . We account for the 
inelastic intermediate states approximating the double virtual Compton scattering 
by the unpolarized forward virtual Compton scattering. The unpolarized proton 
structure functions are used as input for the numerical evaluation of the inelastic 
contribution. Our calculation reproduces the leading terms in the Q 2 expansion 
of the TPE correction and goes beyond this approximation by keeping the full Q 2 
dependence of the proton structure functions. In the range of small momentum 
transfer our result is in good agreement with the empirical TPE fit to existing data. 
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I. INTRODUCTION 

The discrepancy in the ratio of the electric ( Ge ) over magnetic ( Gm ) proton form factor 
(FF) extractions, when comparing the nnpolarized scattering using the Rosenbluth separa¬ 
tion technique [1] with the scattering of a polarized electron beam followed by the recoiling 
proton’s polarization measurement [2-5], or the scattering on a polarized proton target [6], 
spurred a reconsideration of the formalism of elastic electron-proton scattering. The ap¬ 
proximation of the exchange of one photon (OPE) fails at larger momentum transfers, when 
the relative contribution of the processes with exchange of two photons (TPE) increases 
[7, 8]; see e. g. Refs. [9, 10] for reviews. During the last decade, TPE corrections were esti¬ 
mated in different approaches [8, 11-15]. Few direct measurements of the TPE corrections 
were proposed and performed at Jlab [16, 17], at VEPP-3 [18], and at DESY (OLYMPUS 
experiment) [19], with the analysis of the latter currently in progress. 

Besides the importance of TPE corrections at larger momentum transfers, when ex¬ 
tracting the Ge/Gm ratio, recent high precision measurements of elastic electron-proton 
scattering at lower momentum transfer also make it necessary to better quantify the TPE 
effects in order to extract proton FFs. Despite the relatively small absolute correction to 
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the electric and Zemach radii [20], TPE corrections can significantly change the magnetic 
radii [21], and therefore this correction has to be additionally studied. In this context, sev¬ 
eral improved evaluations of TPE box contributions with different intermediate states were 
performed recently for elastic electron-proton scattering [22-27]. 

The importance of the qualitative and quantitative understanding of TPE effects is also 
necessary in view of the 4-7 standard deviations discrepancy in the proton rms charge ra¬ 
dius extraction when comparing muonic hydrogen spectroscopy extractions [28, 29] with 
experiments using electrons [21, 24, 30-32]; see Ref. [33] for a recent review. The TPE cor¬ 
rections constitute the leading hadronic uncertainties in the extraction of the proton charge 
radius from muonic hydrogen spectroscopy. Several recent works [34-39] estimated that 
they account for 10% - 15% of the correction needed to reconcile the proton charge radius 
extractions from muonic hydrogen spectroscopy and electron experiments. 

In the region of very small momentum transfer Q 2 , the TPE corrections can be determined 
model independently. Consequently, each TPE calculation should give the same result or 
at least reproduce some characteristic features in the low-Q 2 limit. The leading term in the 
momentum transfer expansion ( Q 2 ) arises from the scattering cross section of the relativistic 
massless electron on a point charged target in Dirac theory given by the so-called Feshbach 
correction [40]. This result was reproduced later by Brown in Ref. [41] as the leading term 
in the expansion of the TPE correction with the proton intermediate state (elastic TPE). 
Brown also found that the subleading Q 2 In 2 Q 2 term entirely arises from the elastic TPE, 
while the subleading Q 2 In Q 2 term arises from both the elastic and inelastic TPEs. In Ref. 
[41] it was shown that the inelastic TPE correction to the Q 2 In Q 2 term can be expressed in 
terms of an energy integral over the total photoabsorption cross section on a proton target. 
A recent numerical estimate of this term was given in Ref. [42]. 

In this work, we extend the low-Q 2 limit of the inelastic TPE contribution beyond the 
leading Q 2 lnQ 2 term. We express its contribution in terms of the two unpolarized forward 
virtual Compton scattering amplitudes. In detail, the TPE correction is presented as a 
double integral over the virtual photon energy and the virtuality of the unpolarized nucleon 
structure functions. Our results reproduce the known Q 2 lnQ 2 term of Ref. [41] and get 
rid of the hadronic scale A dependence, which arises to the order Q 2 \nQ 2 , thus partially 
accounting for the Q 2 term and higher terms in the momentum transfer expansion. We 
keep the Q 2 dependence in the proton structure functions in order to reproduce the model- 
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independent TPE correction at low momentum transfers, as well as to provide numerical 
estimates in the most accurate way. 

The paper is organized as follows. We introduce the general description of the TPE 
correction to the unpolarized elastic electron-proton scattering in terms of the double virtual 
Compton scattering process in Sec. and reproduce the leading terms in the momentum 
transfer expansion of the elastic TPE contribution in Sec. . We express the leading term of 
the inelastic TPE correction as an integral over the unpolarized proton structure functions 
in Sec. . We show in Sec. that the subtraction function in the forward Compton 
scattering amplitude 7) is negligible for elastic electron-proton scattering. In Sec. we 
describe the details of our evaluation of the inelastic TPE. In Sec. we show that our 
expression reproduces the known result for the leading Q 2 In Q 2 inelastic TPE contribution. 
We present the results of the numerical evaluation beyond the leading terms in Sec. 

Our conclusions and an outlook are given in Sec. 

II. TPE CORRECTION IN TERMS OF VVCS 

The elastic scattering of an electron off a proton target in the OPE approximation is 
described by the helicity amplitude, 

I” 1 = ( k '< h ') (■ k . h ) • * (p', s') r« ( Q 2 ) N (p, s), (1) 

with the electron (proton) initial and final momenta given by k, k' (p, p') respectively, the 
squared momentum transfer Q 2 = — q 2 = — (k — k') 2 , and the unit of electric charge e. The 
kinematics of the process and the electron (proton) initial and final hclicities h, b! (s, s') are 
shown in Fig. . 

The electromagnetic vertex for the transition 7 * (q) + N (p) —>■ N (p + q) is given by the 
Lorentz and gauge invariant form that preserves parity and charge conjugation symmetries: 

n (Q 2 ) = Fo (Q 2 ) r + Fp (Q 2 ) , (2) 

with M the proton mass and Fu and Fp the Dirac and Pauli FFs of the proton, normalized 
to F d (0) = 1 and Fp (0) = k, where k is the anomalous magnetic moment in units of 
e/ (2 M). The Dirac and Pauli FFs are often equivalently expressed in terms of the Sachs 
magnetic and electric FFs as 

Gm — F d + F P , Gp = F d — tF p , 


( 3 ) 
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FIG. 1: Elastic electron-proton scattering in the OPE approximation. 


with r = Q 2 / (AM 2 ). 

The exchange of two photons contributes to the elastic electron-proton scattering through 
the TPE diagram shown in Fig. . 
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FIG. 2: Direct and crossed TPE diagrams. 

The lower blob in Fig. is given by the double virtual Compton scattering (VVCS) 
process on a proton (see Fig. ): 7* (g 1; Ai) + N ( p , s) — » 7* (g 2 , A 2 ) + N ( ps'). The VVCS 
amplitude can be written in terms of the VVCS tensor as 

Tx 2 s’Ms = (qi, Xi) e; (q 2 , A 2 ) • N (p', s') M^N (p, s) , (4) 

where £„, e* denote the virtual photon polarization vectors, N, N the proton spinors, and 
Ai, A 2 (s, s') the photon (proton) hclicities. The initial and final virtual photon momenta q\ 
and q 2 are related to the momentum transfer q and the loop variable q in Fig. by 


~ . q 

qi = q + r 


~ q 

q2 = q ~2- 


( 5 ) 
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For the following, it is also convenient to introduce the averaged lepton (nucleon) four- 
momenta K (P), respectively, as 



FIG. 3: Nonforward VVCS process. 

The TPE correction 5 27 to the unpolarized e~p scattering cross section, due to the graph 
in Fig. in the leading order of the fine structure constant a = e 2 / (47t) can be defined 
as the difference between the measured cross section o and the cross section in the OPE 
approximation cr l7 as 


a ~ <T i 7 (1 + <5 27 ) 


(7) 


The TPE correction <5 27 is given by the interference of the one-photon exchange T l7 and the 
two-photon exchange T 27 amplitudes, 


^ 2 7 — 


2 K [ £ T 27 (T 17 )* 

spin 

E l ^ l7 | 2 

spin 


with the interference term 


2& ( ^T 27 (T 17 )* J = -SR 


47 re 4 f id 4 q 






\spin 


V Lva = TV { | 7 m /t A n Y + Y , A+ A +m O " I ( k + m ) 7 ° ( k’ + m 


(K — q) — m 2 (K + q)~ — m 2 

H >,vo = TV | (p + M) r a (Q 2 ) (p’ + Af) } , 


( 8 ) 


Q2 J ( 2?r P ((<? - D 2 - m 2 ) (G+1) 2 - a 2 ) ’ 


where m denotes the lepton mass, and where /j denotes the small photon mass, that plays 
a role of the IR regulator. 
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Furthemore, the VVCS tensor M ttv reads, in the notation of Refs. [43, 44] 1 , as 

= a Bi {ql ql qi ■ g 2 , q ' P) T? v , 

i£j 

J — {1, 21} \ {5,15,16} , (9) 

where the 18 independent tensors T/ 11 ' were constructed to be gauge invariant, and free 
of kinematical singularities and constraints, following the procedure outlined in Ref. [45]. 
The invariant amplitudes Bi satisfy definite transformation properties with respect to photon 
crossing as well as charge conjugation combined with proton crossing as detailed in Refs. [43, 
44], The lack of knowledge of the amplitudes B\ does not allow one to evaluate the TPE 
correction in the general case. 

In this work, we will consider TPE correction to the unpolarized elastic electron-proton 
scattering in the low -Q 2 limit. The momentum transfer accessed by current electron-proton 
scattering experiments is much larger than the electron mass m. We will consider the 
kinematic range m 2 <C () 2 <C M 2 , ME. 


III. ELASTIC TPE CONTRIBUTION AT SMALL Q 2 


We start our study by checking the low-Q 2 limit for the case of proton intermediate states 
in the lower blob of the TPE graph; see Fig. . In this case, we can compare our result with 
the known result for the elastic contribution. 

The leading terms of the low-Q 2 expansion of the elastic contribution to the TPE correc¬ 
tion 5l° rn in the unpolarized elastic e~p scattering are given by the sum of the IR divergent 
piece 5^, the Feshbach term 8f [40] expanded to the leading order, and the finite logarithmic 
corrections as 


zBorn 

"27 


-A 


$2^ T $F + 


Q 

ME 




+ 0 


Q 2 Q 2 Q 2 \ 

M 2 ’ ME' E 2 ) ’ 


where 5p and are given by 


5 f —t 


omQ 
2 E ’ 



a Q 2 p 2 

nlME^Q 2 ' 


( 10 ) 


( 11 ) 


1 Note, however, that our convention for the photon field four-momenta indices y and v (see Fig. 3) is 
opposite to Refs. [43], [44], which changes the sign of the tensors that are antisymmetric under [i -f-)- v. 









This result was obtained in two equivalent ways in the literature: either by considering the 
low-Q 2 limit of the box graph or by using dispersion relations for the invariant amplitudes 
of the elastic electron-proton scattering; see Refs. [26, 41] for details. 

The correction of Eq. ( ) is given by the graph with two Dirac structures in the 

photon-proton-proton vertices, which are independent of the proton structure. The VVCS 
tensor M pv oint in this case is given by the sum of three nonvanishing terms [43-45]: 


M£ int = a (B p 2 0int T^ + B p ™ nt T p " + B p ™ nt Ttf) 


( 12 ) 


The amplitudes B\™ nt , B p ° mt and the expressions for the corresponding tensors are 

given in Appendix 

The leading elastic contribution to the TPE correction of Eq. ( ) is given by the con¬ 
tribution from the unpolarized VVCS amplitude B^ 0 "' 1 '. Neglecting the subleading terms in 
the momentum transfer expansion (i.e. taking the limit Q 2 <C M 2 , ME , E 2 ), the leading 
terms in the Q 2 expansion for a Dirac point particle can then be obtained as 

C”' J jdhjn+npn+n^n+riQ j (2? + 3 -q(k ■ p) (k ■ q) (p ■ <j) 

-2 (K- q? (p ■ qf - + y) ( m2 ( k ■ i) 2 + Y (p - «> 2 ) 

-^ 2 + t) 2 (k -p )2 + Se (k -^ (p - & j’ (13) 

where we have introduced the propagator notations: 

rip =- \ -, n^. =- —j -, n± =-—,-. (14) 

P (P ± q) 2 — M 2 K (.K±qf-m 2 Q ( q±q/2) 2 -fi 2 

The expression in Eq. ( ) reproduces all terms in the low-Q 2 expansion of Eq. ( ). The 

TPE correction from the spin-dependent structures Bf™ nt T^ and B\™ nt T^ expansion starts 
with higher powers in Q 2 in comparison with the expansion of the unpolarized structure 
B'po'mirppv contribution to TPE, and consequently the leading Q 2 term in the expansion 
coming from the graph with two vertices is reproduced in a correct way with Eq. ( ). 

As an example, we describe the derivation of the Feshbach correction [the second term in 
Eq. ( )] in Appendix 
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IV. INELASTIC TPE CORRECTION IN TERMS OF THE FORWARD 
VIRTUAL COMPTON SCATTERING AMPLITUDES 


The VVCS amplitudes B\, B 2 , B%, B$, Big contribute to the unpolarized forward virtual 
Compton scattering tensor. 

The elastic intermediate state contribution from the unpolarized VVCS amplitudes is 
described by two nonvanishing amplitudes, 


jz^Born _ 


V 


M hi ■ ®0 2 - (2 Mv) 2 

(91 - 92 ) 


{F P (Ql) F P (Ql) 


(F d (Qf) F P (Qf) + F P (Ql) F d (Ql) + F P (Ql) F P (Ql)) j , 

(15) 

( Fd m Fd (QI) ~ ^ Fp(Ql) Ff 
with Ql = —qf, Ql = —ql, the kinematic relation (gi • q 2 ) — — (ff 2 — and notations 


B. 


Born 


AK = (P-q), Q 2 = 


(17) 


These amplitudes contribute to the VVCS tensor of Eq. ( ) with the following tensor 
structures: 


rr = - (31 ■ 32) u + 3i 321 (!8) 

= -4 (P ■ qf g> w - 4 (qi ■ q 2 ) P»P V + 4 (P ■ q) {P v q£ + P^qf) . (19) 


The other unpolarized amplitudes vanish in the Born approximation, i.e. 

Bf orn = 0, Bf orn = 0, B*° rn = 0. (20) 

When evaluating the TPE correction, one can simplify the calculation by using explicitly 
gauge invariance, i.e. q 2 L^ wa = 0 and q\L liua = 0 for the virtual photon momentas q\, q 2 . 
Consequently, we can use the identities: 


qiL^ a = q^L 


fivon 


0.2 L iiv& Q L 


fiua • 


( 21 ) 


Approximating the arguments of proton form factors as: Ql ~ Q\ ~ Qf — Q 2 / 4, the Bf orn 
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and B 2 orn contributions can be obtained from the effective near-forward VVCS tensor, 



( 22 ) 


with the Born contributions to forward VVCS amplitudes T® orn and T^ 0 ™ given by 


Tf orn (z/,Q 2 ) = 



Q a G 2 m ( Q 2 ) 





(23) 


(Q 2 ) 2 - 4M 2 i/ 2 


T 2 ° rn (c Q 2 ) = 4 MQ 



( 1 Q 2 ) 2 - 4M 2 z/ 2 


(24) 


Note that in Eq. ( ) the second argument of the VVCS amplitudes is given by 

replacing in Eqs. (23) and (24) Q 2 —» Q 2 — Q 2 / 4. 

The expression of Eq. ( ) gives the correct forward limit, when q — 0. For a Dirac point 

particle, we checked that Eq. ( ) corresponds with the exact result for the unpolarized 

VVCS tensor when replacing F n ( Q 2 ) —* 1 and F P ( Q 2 ) —>■ 0 in Eqs. ( ) and ( ). At low 

Q 2 , it leads to Eq. ( ) and gives the leading terms in the low-Q 2 expansion of the elastic 

contribution to <5 27 of Eq. ( ). We call the approximation of Eq. ( ) the near-forward 

approximation. 

We now turn to the inelastic contributions to the unpolarized forward Compton ampli¬ 
tudes Tx, T 2 . They are expressed in terms of the VVCS invariant amplitudes Bi in forward 
kinematics by 


(s, Qf = a (<§ 


2 B! - 4M-V-B, + Q'Bi - 4MvQ 2 B t ) , 


(25) 

(26) 



In this work, we choose the tensor form of Eq. ( ) to describe the inelastic TPE correction 

and only keep the amplitudes B 1 and B 2 in the region of small momentum transfer. 

We describe the VVCS invariant amplitudes B, (g 2 , q 2 , qi ■ q 2l q ■ P ) as 


Bi (g 2 , q 2 , gi-g 2 , g-P) ^ B t 



where we use the approximation g 2 « q 2 pe —Q 2 + Q 2 /4 for the VVCS amplitudes. 
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The near-forward approximation allows one to obtain the first two terms of the inelastic 
TPE expansion coming from the proton structure functions F\ and F 2 : 


zinel 

°2 1 


a(E ) Q 2 lnQ 2 + b(E ) Q 2 . 


(28) 


Nevertheless, we keep the Q 2 dependence in all kinematic factors, but we do not pretend on 
the validity of our approximations beyond the expansion of Eq. ( ) due to the contribution 

of double virtual Compton amplitudes besides Bi and B 2 . The near-forward approximation 
of Eq. ( ) is valid only in the region of small momentum transfer Q 2 . In Sec. we will 

explicitly study the range in Q 2 over which such expansion provides a good approximation. 

We obtain the TPE correction substituting the near-forward approximation of the VVCS 
tensor of Eq. ( ) into the general expression for the cross section correction 5 27 of Eq. ( ): 


fiinel _ 4 _(1_ £ )Ge 

27 - sG 2 + rG 2 M* 


X 


T E - t - I 

idYq 


(27p 4 Q Q 

Q 2 


U Q U Q X ■ P) ((K ■ py - ^-P 2 


Q 2 „,\ n 


K 


n 


K r 


M 2 


Q 2 -¥ 

/ 


((P ■ q) (K -q)P 2 -(P- qf (K ■ P)) Uk+ ^ k T, 


M 2 


- m2 Q 2 
V,Q -T- 


Q 5 


- —r 


j->2 , {P-qfQ 2 

Q 2 

4 


\ 


V 

—2 (P ■ q) 


7 


(k ■ p) (. k ■ a) - ^ (p ■ a) 


K r 


- m2 Q 2 

”,Q ~ ~r 


Q 2 - 2 ? 

q 2 -¥ 


(K ■ PY - —P 2 
y ' 4 


Q 2 ^ - n 


M 2 


K 


~ m 2 ^ 
v,Q -~r 


.Q 2 + % 


((A- ■ P) (A. a - ^ (P ■ a) (n* - n+) t, U Q‘ 


Q 2 -f 

with the photon polarization parameter: 

16 (K ■ P) 2 - Q 2 ( Q 2 + AM 2 ) 


91 

4 


(29) 


£ = 


(30) 


16 (A • P) 2 + Q 2 (Q 2 + 4M 2 )' 

The imaginary parts of the forward virtual Compton amplitudes T 1; T 2 , by means of the 
optical theorem, are directly related to the unpolarized proton structure functions F\ and 
F 2 by 


= 4M Fl ' 


3T, = — P 2 . 
Au 


(31) 


To obtain the real parts of the forward Compton amplitudes, we express the dispersion 
relations (DRs) in terms of the invariant mass of the intermediate state variable W 2 in 
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order for the boundaries of the integration regions not to depend on Q 2 . This allows us to 
easily interchange the W 2 integration with the q integration. The DRs for forward virtual 
Compton amplitudes are given by 


T,(S,Q 2 -^ 


= Ti[0,Q 2 —^ 


2 

71 


e 2 Mv 2 F x (w 2 ,Q 2 - dW 2 


w 2 , 

vv thr 


v,Q ~ ~r ) = ~ 


(w 2 - P 2 + Q 2 ) ((P + qf - W 2 + ie) (( P - qf - W 2 + ie) ’ 
e 2 MF 2 (w 2 , Q 2 - Of) dW 2 


vr J ((P + q) 2 - W 2 + ie) ((P - qf - W 2 + ie) ’ 


(32) 


w\ 

vv thr 


with the pion-proton inelastic threshold, W 2 hr = (M + m 7r ) 2 ~ 1.15 GeV 2 , where rn, denotes 
the pion mass. The unsubtracted DR for the Ti amplitude is divergent due to the Regge 
behavior of F\ structure function. Consequently, we use a once-subtracted dispersion relation 
for this amplitude with the conventional subtraction point z> = 0. 

We find the contribution from the subtraction function T| ^0 ,Q 2 — to be negligible 
in electron-proton scattering experiments; see the following Sec. for details. The contri¬ 
butions from the unpolarized proton structure functions F\ and F 2 to the TPE correction, 
which we denote by 5^ and 5^, respectively, are given by 


°27 


fi F 2 

°27 



w 2 , 

vv thr 



W 2 , 

vv thr 


mg n+n^B (n K - n+.) (p ■ gf f, (\v\ q 2 - g) 

(2?r) 4 f w 2 _ p2 + Q2\ ((P + ,J|2 _ H'2 + j (,p _ _ U™ + fc) ’ 

(33) 

id“« n O n Q (<? ( n A' + n t) + D ( n * - n b) F ^ ( W 't Q 2 ~ ¥) 

----- — (34) 

(2vr) 4 ((P + g) 2 - W 2 + is) ((P - qf - W 2 + ie) 
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with 

A 

C 


b = ShhL (q2 (p . a - 4 (x. p) . a), 

V 4 

Q 2 


8e 2 (l-e)G E 
M 2 + tG 2 m ’ 

2 (if ■ p) ((a- ■ p) 2 - 5!p 2 ) + (p ■ a «A- ■ a -p 2 - (p ■« (a ■ p )), 

N ' Q 2 


D = P z + 


2l (p-q) 2 Q 2 


Q 2 ~ Q r 


P) / 


(P-P) (X-g)-^(P-g) 


/02 _ 3Q 2 / 

- 2 (A-a|^((A-p) 2 -^p 2 ]. 




(35) 


V. T, SUBTRACTION FUNCTION CONTRIBUTION 


Before discussing the inelastic contribution to <5 27 , we first discuss the subtraction term in 
the T x amplitude regardless of the near-forward approximation of Eq. ( ). The subtraction 

function Tx (0, Q 2 ) is expressed in terms of only two VVCS amplitudes Pi and P 3 as 

Tx (0, Q 2 ) = a (Pi (Q 2 ) Q 2 + P 3 (Q 2 ) Q 4 ) , (36) 

with the relevant VVCS tensor structures given by 


rjiflV 


rj-}flV 


- (<Zi • <? 2 ) + «, 

qlqlg^ + (?i • p) <?r<?2 


9l + <? 2 2 


( 91«1 + Q2Q2) + 


9? 




(<?r<?r 


9^2 ) 


The TPE correction due to the subtraction term, arising from the amplitudes Pi, 
is obtained from Eq. ( ) as 


(37) 

(38) 

B 3 , 


rsubt 

h'y 


id 4 q 

Q 2 


Pxn+n- ( g 2 + ^ ) + p 3 ) n+n 


Q Q 


12 , Q' 


L K LL K 


X (P ■ g) (tf ■ q) - (K ■ P) (K ■ q) z j = h + / 2 - 
We express the loop integral of the first term of Eq. ( ), denoted by It, as 


(39) 


cy f iif-q 
1 4 J (2n) 4 

Q 2 f ii l q 

16 J (2jt) 4 


BinAK ( q 2 + % ) + b 3 ) (p ■ a (a ■ a iyn 


1 K 


Q- 


Bin+n 0 e 2 + ly + S 3 (P ■ a (n* - n+ ) = a K (k ■ p) , (40) 
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where we have defined ax as 

° kA ’" = m / 0 ( BlE S n 5 (« 2 + t) + B3 ) ®" ( n « “ n t) ' < 41 ) 

On the other hand, since K 2 = Q 2 /4, we can rewrite the integral I\ as 

h = (K (siU+iq (<5“ + ^) + fla) (A- ■ n+n*, (42) 

which exactly cancels the integral I 2 from the second term in Eq. ( ). Consequently, the 

subtraction function contribution to unpolarized elastic electron-proton scattering vanishes 
in the limit of massless electrons. 

It is instructive to study the subtraction function contribution to the TPE amplitude T 2 " 7 
of elastic lepton-proton scattering, accounting for a finite lepton mass. The amplitude is 
given by 


T 2 7 = _ 2 / L^N(^, s ')M^N(p, s ) 

' J C 2 ^) 3 ((g - |) 2 -/x 2 ) ((?+|) 2 -m 2 )’ 


(43) 


with the leptonic tensor 


tuv ~ni u\ I u K-q + m u v K + q + m . 

= u (k, h!) Y - Y - 7 + 7 --U-7 u ■ 

V n (K - qf -m 2 (. K + q ) 2 - m 2 1 V J 


We study the TPE contribution due to the VVCS amplitude Bi first. The — (gi ■ q 2 ) g^ 
term contribution is given by 


zd 4 g 




Bi ( NN ) • 2m (q 2 + ^ ) uu + 4 (K • q) hg~« g~ 2 + ^ 11+11^11+1^. (44) 


Q 2 


Using the gauge symmetry in the tensor T(' M the q±q 2 term contribution is given by 

id 4 q 


e 4 Q 2 


Denoting 


e 4 Q 2 


(2tt) 


id 4 g 


B, (JVJV) • 2 (K ■ q)uqu) n+I^IIJiq. 


(45) 


Bi ■ ( 2 (A- ■ q)r) n+n^n+n^ = 


(46) 


the contribution of the gfg^ term is given by 


ma'x (NN~) ■ (uu). 


( 47 ) 
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Note that, according to the symmetry properties of the VVCS amplitudes [43], the 
expansion of Eq. ( ) also contains no q M term for the case of VVCS amplitude 
Pi (f/i, ?|, (q\ ■ Q'i ), (P • q)) subtracted at an arbitrary point (P ■ q) = Mi/ 0 . The ampli¬ 
tude P 3 contributes through only one tensor structure qfq^g^ in the similar way as the 
amplitude Pi contributes through the structure — (q± ■ g 2 ) . 

We conclude that the Ti subtraction function contributes to the (uu) ■ (NN) term in the 
elastic lepton-proton scattering amplitude. This contribution is suppressed by the lepton 
mass 2 . In the language of effective held theories, the chiral (or axial) symmetry on the 
lepton side forbids the (uu) ■ (NN) type structure for massless electrons, and therefore the 
contribution of the subtraction function vanishes in the massless lepton limit. 

VI. EVALUATION OF THE INELASTIC TPE CONTRIBUTION 

In this section, we describe the strategy to evaluate the inelastic TPE correction S^ el , 
which we will express as a two-dimensional integral over the unpolarized proton structure 
functions F\ and P 2 ; 

= J <w 2 d< 3 2 [»', fir 2 q‘) /■, (u J . (j- - Tj 

+w 2 (w\ q 2 ) f 2 (\v\ & - ^P) |. 

(48) 

To obtain the result of Eq. ( ), starting from Eqs. ( ) and ( ), we will perform the 

integration, using two choices of Breit frames with the aim to cross check the method and 
its application. In a first choice, which we denote as P-frame, the kinematics of the external 
particles is defined by 

K= (iP o ,0,0,|i?|), P = P( 1,0,0,0), g = Q (0,1,0,0). (49) 

In a second choice, which we denote as it-frame, the kinematics of the external particles is 
defined by 

K = K (1 ,0,0,0), P= (P 0 ,0,0,|P|), q = Q( 0,1,0,0). (50) 

2 That is equivalent to a subtraction in the dispersion relation for the E 4 invariant amplitude, in the notation 
of Ref. [46] . 
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In either of these frames, we evaluate the weighting functions w i, W 2 by performing 
a Wick rotation in Eqs. ( ) and ( ) in the complex q 0 variable plane. The resulting 

TPE correction is given by the sum of the integral along the imaginary q 0 axis and the 
contributions from the poles which are crossed by the integration contour. In addition to 
the photon and lepton propagator poles, there exist the hadronic poles coming from the 
dispersion relation propagators: 

H ( q±P) 2 -W 2 

The pole positions are shown in Fig. . 




FIG. 4: The position of the % poles for different propagators. The lepton poles (11//) and hadronic 
poles (II//) contribute only in a limited region of integration variables. 

The integration contour does not cross the photon poles 11^. The expressions of Eqs. 
( ) and ( ) are symmetric with respect to the change of integration variable q —» —q. 

Exploiting this symmetry, we only need to calculate the residues of the upper half plane 
poles and double the result. The leptonic pole 11^ and the hadronic pole 11^ are moving 
poles; they contribute in the limited region of W 2 , Q 2 variables. We show the corresponding 
regions in Fig. 5 for the limit of low-Q 2 scattering. 
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0 


P(l, 0, 0,0) frame 

4E 2 


K(l, 0, 0, 0) frame 

_Q! n of 

4 U 4 

mmmdUWM . 


► 



16M 2 E 2 

Q 2 


► 



integral along the imaginary axis 
K pole contribution 
P pole contribution 


FIG. 5: Integration ranges for the different contributions to ■ 

We provide the expressions for the regions of these poles contributions and weighting 
functions Wi,w 2 in Appendix . 

VII. LEADING Q 2 lnQ 2 INELASTIC TPE CONTRIBUTION 

In this section, we describe the way to obtain the leading Q 2 In Q 2 inelastic TPE contri¬ 
bution in terms of the total photo absorption cross section a? [41] in our approach. For the 
region of small momentum transfer, we exploit the approximation [47]: 

Fi (W 2 , Q 2 ) « {W 2 ) , (52) 

F 2 (w 2 ,Q 2> ) « ^a T {W 2 ). (53) 

In the P-frame, the leading correction is coming from the leptonic pole, whereas in 
the /l-frame, it originates from the hadronic pole. We present its derivation in the K- 
frame. The pole q 0 = P 0 — ^ J ^P ~ <tj ~ IP 2 contributes in the invariant mass region 
W 2 < Pq ~ 4 (K ■ P) 2 /Q 2 . For the leading contribution, we consider the W 2 integration 
up to infinity. For the case of small momentum transfer values, the integration region in the 
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Q 2 variable is given by 


~ 2 16 {K-Pf 


a 0 Q z <Q 2 < 


Q 2 


with 


ot o — 


W 2 - M 


2 \ 2 


(54) 


(55) 


4 (K ■ P) 

We can replace the upper integration limit by infinity, but the approximation of Eqs. ( ) 

and ( ) in terms of the photoabsorption cross section is valid only up to some hadronic scale 

A, that reproduces the Q 2 behavior of unpolarized structure functions 3 . In the present work, 
we will go beyond such approximation and directly use the unpolarized proton structure 
functions F\ , P 2 with their Q 2 dependence as an input. The integration region Q 2 > A 2 
will not contribute to the term proportional to Q 2 lnQ 2 , and therefore we integrate up to 
the squared hadronic scale A 2 . 

The logarithmic term comes from the region Q 2 Q 2 . Accounting also for the pole 

condition (P ■ q) = (p 2 — Q 2 — W 2 ) /2 » Q 2 , we get 


rFl 

0 Orv - 


2 Q 2 e 2 


27 (K ■ P) 


/ dW 2 


d 3 g Kq 0 (n~-n+) F i{w 2 ,Q 2 


(2vr y 


Q 4 


Po-q o 


(56) 


w 2 , 

yv thr 


C 2 = 4Q 2 e 2 9f? / dW 2 


i 3 q (p ■ ?) (n« - n+) - (k ■ p) (n* + n+) ^ [w 2 ,Q 2 


W? hr 


(2tt) 


Q 4 (W 2 - P 2 + Q 2 


Po ~ Qo 


(57) 


With account of the proton structure functions approximation of Eqs. ( ) and ( ) the 

leading term of the TPE correction expansion is expressed as 

OO A 2 

Si 1 = - Q4 [ ^a T (W 2 ) [ dQ 2 ^-^ (Li-Lt), (58) 


"2 7 


W2 

— 


(.K ■ Pf 
Q 2 K 


167T 3 


w 2 

yv thr 


OtoQ 2 

A 2 


Q 4 


27 (K ■ P ) 


3? 


m F Ti w 2 ) A Q 2 Aiim±S^FMF^il m 


47T 3 


w 2 

tv thr 


«0 Q 2 


Q 2 (P ■ q) 


3 Note that the numerical evaluation of Ref. [42] corresponds with a specific choice of hadronic scale, 
A 2 = (W 2 - M 2 ) 2 / (M 2 + 2 ME). 
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with the following integrals: 


n M 

% 


% 


Lt = 


dg 0 n| = 


dg 0 _ 1 W 2 - P 2 =f 2 (K • P) 

Q 2 TQ%~ ± Q ln w 2 -P 2 


(60) 


„-M 


=M 

?0 


= 


g 0 dg 0 n^ = 


9odg 0 

< 2 2 =F < 5 <?o Q 


Q 2 ( W 2 — P 2 =)= 2 (A' • P) 2 (AT ■ P) 
^ — In-——--- 1 ± —-- 


W 2 - P 2 


W 2 - P 2 


% 


M 


2 (k ■ p ) g 2 


W 2 - P 2 Q 

Performing the Q 2 integration, we obtain the leading logarithmic contributions, 


(61) 


cfi _ (91 

^ 87r 3 n V A 2 


Q\ (Q 2 


d W 2 (W 2 -P 2 1 2(K ■ P) + P 2 -W 2 , , o \ , . 

1 In ^^ + 1 a T (W 2 ), (62) 


w 2 , 

vy thr 


ME V 4 ME 2 (A' • P) - P 2 + IP 2 


A 2 


Q\ (91 

4vr 3 1 1 A 2 


dlP 2 


■ 2\2 


In 


(IP 2 - P 2 ) 




ME (2 (AT • P) + P 2 - IP 2 ) (2 (A' ■ P) - P 2 + IP 2 ) 


dlP 2 , 2(K ■ P) + P 2 -W 2 /TI72n 
I n ^— ^a T (W 2 ), 


a T {W 2 ) 


w? hr 


W 2 -P 2 2 (AT • P)-P 2 + W 2 


(63) 


which sum up to the known expression of Ref. [41] . 

When also accounting for the Q 2 terms in the expansion of Eq. ( ) for the TPE correc¬ 
tion, the hadronic scale A dependence in Eqs. ( ) and ( ) drops out. 


VIII. RESULTS AND DISCUSSION 

Having verified the leading Q 2 In Q 2 contribution for the inelastic TPE, we next discuss 
the numerical evaluation of Eq. ( ) including the Q 2 dependence of the structure functions. 

As a numerical check, we evaluate the weighting functions w\ and W 2 both in the K- 
frarne and P-frame, in Figs. . When adding the integral along the imaginary axis with 
the pole contributions, we checked that we obtain the same result in both frames. Despite 
the finite Q 2 ranges of the poles contribution, which are distinct in the K- and P-frames, 
the resulting weighting functions are continuous. The boundaries of the A"-pole region and 
the lower boundary of the P-pole region, shown in Fig. , are covered by Figs. . The 
weighting functions W\ and W 2 have a singularity at Q 2 = Q 2 / 4, when the photons are on 
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their mass shell in our approximation, i. e. qf = q 2 = 0. The weighting function W 2 has a 
discontinuity in the first derivative at Q 2 = 0, where the integral along the imaginary axis 
starts to contribute. 



FIG. 6: The weighting functions of Eq. ( ), with the absolute value plotted and where the sign 

labels show the corresponding signs of w\ and W 2 for e~p scattering. As a check, the weighting 
functions are evaluated in two frames, yielding exactly the same result. The kinematics is chosen 
as shown on the figure. In the upper panels, the weighting functions w\ and W 2 are shown in the 
region — < Q 2 < Hr and in the lower panels for Q- < Q 2 < 8 E 2 . 

We perform the Q 2 integration hrst and express the resulting TPE correction S^ el in 
terms of the W 2 integral as 

f(W)dW 2 . (64) 

In Fig. , we present the result for the W 2 integrand for different inputs of proton structure 
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functions and compare our full calculation with the approximation of Eqs. ( ) and ( ) 

with A « 0.6 GeV. To describe the Q 2 dependence of the unpolarized structure functions, 
we use the empirical fit performed by Christy and Bosted (BC) [48] , while for the logarithmic 
approximation in the Fig. 7, we also use the SAID Partial-Wave Analysis Facility [49] for 
the photoabsorption cross section. The BC fit is valid in the region 0 < Q 2 < 8 GeV 2 , 
M + m,r < IF < 3.1 GeV. This fit is not very accurate in the threshold region near 
Q 2 = 0, although it is still compatible with the error bars of the photoproduction cross 
sections. We checked that difference between using the SAID and BC fit in the region 
IF < 1.15 GeV and 0.9 x Q 2 /4 < Q 2 < 1.1 x Q 2 / 4 only amounts to the relative change in 
S^ el at the 3%-4% level for the kinematics shown in Fig. . For the result of 5%F l beyond 
the Q 2 lnQ 2 approximation, we use the BC fit for the region 0 < Q 2 < 12 GeV 2 , where 
the integrand behaves in a smooth way. The relative contribution to 8 l E el from the region 
8 < Q 2 < 12 GeV 2 is smaller than 0.1%. We perform the IF 2 integration up to W = 4 GeV. 
The extrapolation from IF = 3.1 GeV (upper bound of the BC fit) to IF = 4 GeV leads 
to an additional relative contribution to 8 l £f l of less than 2%. We have also checked on 
the SAID parametrization that the region IF > 4 GeV has a relative contribution to 8l£ el 
of less than 2%-3%, when interpolating the SAID parametrization to the Regge behavior. 
The main inelastic TPE contribution is given by the nN -channel. The singular peak at 
IF 2 = M 2 + 2 ME corresponds to the quasireal photon singularity (when both photons in 
the two-photon box are quasireal and collinear with either lepton). This singularity appears 
only for the beam energies above the pion threshold. 

In order to clarify the validity of the inelastic TPE estimates, we study numerically the 
low-Q 2 expansion of the 5™ f el coming from the F\ and F 2 structure functions. In Fig. 8, we 
present the ratio between the TPE correction S™ f el ( Q 2 ) and the low-Q 2 fit S™ el >f lt ( Q 2 ) in the 
form of Eq. ( ) for the energies of available data. We compare the Q 2 /E 2 and Q 2 / ( ME) 

expansions. We perform the fit in either the variable Q 2 /E 2 or Q 2 / ( ME ) in a range which is 
100 times smaller than the range displayed in Fig. . The comparison of our full calculation 
with such a fit over an extended range provides us then with a quantitative argument on the 
Q 2 range where such an expansion holds. If we use as a criterion that the full calculation stays 
within 10% of the fit, we can see from Fig. that for energies corresponding with available 
data an expansion of the type of Eq. ( ) holds for Q 2 < E 2 and requires Q 2 < ( ME ) /5. 

We expect the same type of expansion for the TPE contributions from other amplitudes. 
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FIG. 7: W 2 integrand / of Eq. ( ) for two different external kinematics. The integrand with 

the unpolarized proton structure functions from BC [48] is shown by the blue solid curve. The 
leading logarithmic correction in the approximation of Eqs. ( ) and ( ) with A = 0.6 GeV and 

the photoabsorption cross section from the fit of BC (SAID [49]) is shown by the blue dashed (the 
red dash-dotted) curves. The dominant 7rA r -channel contribution is shown for the SAID fit by the 
red dotted curve. 

The resulting inelastic TPE correction as a function of Q 2 for the beam energy E = 
0.18 GeV is shown in Fig. . We compare the Feshbach term for pointlike particles, the 
elastic TPE correction in the model with dipole form factors [26], the approximation of Eqs. 
( ) and ( ) with A m 0.6 GeV and the result of Ref. [42], The difference between both 

Q 2 In Q 2 curves comes from the term of order Q 2 and appears due to the different choices of 
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FIG. 8: Q 2 dependence of ratio of the inelastic TPE correction to the low -Q 2 fit of the form 
6'. ineijit __ a q2 i n Q 2 _j_ 5 (E) Q 2 ' The ratio is shown as function of Q 2 / (ME) on the left panel 
and Q 2 /E 2 on the right panel. 

A. The approximation of Eqs. ( ) and ( ) implies the same hadronic scale for all inter¬ 

mediate states, while the hadronic scale in Ref. [42] depends on the intermediate state. We 
present our results for the total TPE correction including the Q 2 dependence in the structure 
functions and extrapolate them to the region Q 2 > ME. Increasing the momentum transfer 
the inelastic TPE correction shows a clear departure from the Q 2 In Q 2 term, reducing the 
latter value. With account of the inelastic intermediate states the TPE correction in the 
low-Q 2 region comes closer to the Feshbach correction in comparison with the elastic TPE 
correction only. The inelastic TPE correction has the same order of magnitude and the 
opposite sign in comparison with the proton form factor effects in the elastic TPE. 

We compare the TPE corrections as a function of the e variable for Q 2 = 0.05 GeV 2 
(Q 2 = 0.25 GeV 2 ) in Fig. ( ). The inelastic proton excitations compensate the proton 

form factor effects, and the resulting TPE correction comes closer to the Feshbach term. For 
the small momentum transfer Q 2 = 0.05 GeV 2 , where we expect the validity of the near¬ 
forward approximation for e > 0.7 or Q 2 / (ME) < 1/5, our calculation is in good agreement 
with the empirical TPE fit of Ref. [21] in the region £ > 0.35 — 0.4 as one notices from 
Fig. . One should definitely account for contributions beyond two unpolarized proton 
structure functions for smaller £ values. We show the region of small £ with the aim to 
illustrate the characteristic features of our calculation. Increasing the momentum transfer, 
as shown on Fig. 1, the predicted TPE correction is found to be in reasonable agreement 
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with the empirical fit of Ref. [21], confirming the proton charge radius values extracted with 
this TPE correction. 


IX. CONCLUSIONS AND OUTLOOK 

In this work, we have evaluated the leading terms in the momentum transfer expansion 
due to the inelastic TPE contribution to the unpolarized electron-proton scattering. We 
expressed this inelastic TPE contribution as an integral over the VVCS amplitudes. The 
VVCS tensor was then approximated by the two unpolarized forward Compton scattering 
amplitudes. The latter are obtained as a dispersion integral over the unpolarized proton 
structure functions. This approximation allows one to reproduce the existent results for the 
low momentum transfer expansion of the TPE correction and goes beyond the Q 2 In Q 2 term 
in such an expansion. Our results for the kinematics of the MAMI experiments show that 
the inelastic contribution has the same sign as the Feshbach correction. With account of the 
inelastic contribution, the TPE correction comes closer to the Feshbach term for pointlike 
particles in comparison with the elastic contribution only. In the limit of low Q 2 , the 
TPE correction is in good agreement with the empirical TPE fit of the Al Collaboration at 
MAMI [21] at small momentum transfers. This agreement confirms the value of the extracted 
electric charge radius from the data of Ref. [21]. Increasing the momentum transfer, the 
TPE correction starts to deviate from the data fit, still following its shape. To further test 
the region of applicability of the near-forward approximation, a detailed dispersion relation 
analysis of separate channels contributions (e.g. the 7rV-charmel) for nonforward kinematics 
will be needed. Such a comparative study will be performed in a forthcoming work. 


Appendix A: Some VVCS amplitudes and tensor structures 

The three nonvanishing VVCS amplitudes for the case of two Dirac couplings in the 
photon-proton-proton vertices in Fig. are given by 


Bf nt = -^iy, 


E 


point 

10 


B 


point 

17 


2 M 

(P ■ q) 


M 


IIpIIp, 
IIpIIp, 


(Al) 
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where the relevant VVCS tensor structures have the following form: 

Tr = -4 (P ■ q) 2 cr - 4 (gi ■ q 2 ) P»P V + 4 (P • q) (P"gf + P"g£), 

= -8 (g, • g 2 ) P^P" + 4 (P • g~) (PV + P^g 2 ") + 4M (q, ■ q 2 ) (P^ + P u Y) 

-AM (P ■ q) (q?Y + fi-f) - 2 (P ■ q) {(g 2 V ~ 9? Y) 7-9 7-9 (tfV* - 9^)} 

-2 (g x • g 2 ) (-P ' 9) (t'V* - YY) + M (9i ' 92 ) {(t'V 1 - 7^7") 7-9 + 7-9 (t'V - YY)} , 

V = -4 (P • g~) + 2 (PV + P^g 2 ) + 4M<T 7 .g - 2M (gV + g 2 V) 

+ W - 9i7") 7-9 - 7-9 (tfV - 9i7")} + (9i ■ 92 ) (tV* - 7V) • ( A2 ) 


Appendix B: Evaluation of the forward limit of the Feshbach correction 


The leading term in the low momentum transfer (Q) expansion for S 2j comes from the 
following term in Eq. ( ): 


V -> 


7h£0!f! J |JnMn TO ^ + f 


E 


(Bl) 


We perform the integration in Euclidean space in this Appendix. The Euclidean coordi¬ 
nates q E can be expressed through the Minkowski coordinates (u, q] by 

(B2) 


qE q E , qE) =qE\ -iv, g , 


{a E ■ b E ) = - (a • b ), 


where we use the index ”E” for the notation of vectors in Euclidean space. 

The leading term of the Feshbach correction can then be cast into the form 


0 MQ 2 e 2 f Yq E 


°27 “ 


( n p + n+) (nj + iq) n Q , 


(B3) 


E J (2 n) 

with ITq — 1/ p|. + 7 + fi 2 ^j . We can expand Ilp^ as Gegenbauer polynomials [C n (a)] 
generating functions, 
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(B4) 
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with 
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Z » = 
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Q 


4 Mx 


l - x + J (l- x 2 ) + 


2 , 16 M 2 

-x A 


Q 2 


1 + X 2 + fl 2 — yj (1 + X 2 + fi 2 ) 2 — 4x 2 
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2x 

y 2 = 4 y 2 /Q 2 . 


(B5) 


Using the Gegenbauer polynomials value C n (0) = (—i)"/ 2 . (i -|_ (—l) n ) /2 and the orthog¬ 
onality relation for vectors q, x , y in Euclidean space 


f / . 2v t 2 

/ dff C n C m ( qy ) ■ ~7$m,nCn \xy) i 

J n + 1 

the integral of Eq. ( ) simplifies to the following expression: 
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(B7) 


For x —y 0, the product ZmZ^ 4 ^ 1 — is IR finite. Numerically, the result of the 

integral does not depend on small jj and y 2 values, so we can neglect terms of order y 2 , Q. 
The resulting TPE correction is given by 


r0 2a Q I x ( 1 + x 
tC = —— / da;--- In 


27 71 E 


1 + X 2 


+ O (Q 2 In Q 2 ) = ^ + O (Q 2 In Q 2 ) , (B8) 


reproducing the low-Q 2 behavior of the Feshbach correction; see Eq. ( ). 


Appendix C: TPE correction in terms of the unpolarized proton structure functions 


In this Appendix, we present the expressions for the pole contributions as well as the 
contribution from the integral along the imaginary axis to the weighting functions W\, w 2 , 
which appear in the TPE correction of Eq. ( ). 

We first present the results in the P- frame, defined by Eq. ( ). The contribution to the 

weighting functions Wi, w 2 arising from the leptonic pole % = Ii 0 — \ (k — q') + m 2 in 


















27 


the P-frame is given by 
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with coefficients P and C + D following from Eqs. ( > 5 ) as 
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In Eq. ( ), Qo stands for the sum of two integrals with either ± signs: 
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The integration regions are given by 
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We do not consider the hadronic pole in the P-frame. It contributes only in the region 
of large momentum transfer: 


Q 2 > 8 Mmv (l + |^) « 1.09 GeV 2 . 


(C6) 




















The contribution to the weighting functions arising from the integral along the imaginary 
axis in the P-frame is given by 
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with the notations 
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with 


c = Q' 2 -Q- + 2K 0 q 0 , g 2 = AK 2 Q 2 sin 2 0. 


(CIO) 


The integral along the imaginary axis contributes in the range 0 < Q 2 < oo. 

As a check of our calculations, we also provide the expressions evaluated in the K -frame, 
defined by Eq. ( ). The leptonic pole go = K — s/qP^rm 2 contribution to W\, w 2 evaluated 

in the A-frame is given by 
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with the notations 
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in terms of the integrals of Eqs. ( ), where 
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The contribution to w i, vj 2 arising from the hadronic pole go 
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a P 2 _g^_1 C + P 1 

2vrcG| + tG 2 m M 2 (K ■ P) 2 + M 4 r (1 + r) |P| Q 2 + ^ (P ■ g) A 
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with the notations 

b = ( (ir ' p) K P - L t) - t (F ' ® ( L » - L « + )) • 

c = (K ■ P) (2 (K .pf-9Pl- (P . (X- + L +) 

+ ^Si (F '® p2 hr + ^), 

V 4 


D 


P 2 + 


(P-g) 2 g 2 


2 _ Of 

4 


(K . P) A- (Lf - Li) 


Q 2 --£ 


■ 2<P '«|?T^h^'-P) 2 


Q 2 P 2 


(Lo~L o + )) 

it)- 


and master integrals 



(C16) 


(C1T) 


The integration region is given by 




< Q 2 < (lei + v^yhv 

< bb 2 < P 2 . 


(C18) 


The contribution to the weighting functions arising from the integral along the imaginary 
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axis in the A'-frame is given by 

oi G e 


WilW^Q 2 ) = 


P 2 


Q 4 


Q- 


87r 2 eG% + tG 2 m M 2 (.K • P) 2 + M 4 r (1 + r) Q 2 - ^ 

7T 

A sin 2 


x3« 


Q 2 — a + 2iKQ cos -0 


W2 (^ 2 ,g 2 ) = 


a G r P 2 


g 4 


g 2 


8tt 2 eG\ + rGlf M 2 (AT • P) 2 + M 4 r (1 + r) g 2 + ^ 

f (c + D) sin 2 0d0 
xK 


g 2 - ^ + 27A'g CO s0 


(C19) 


with the notations of Eqs. ( ), the integrals of Eqs. ( ), and notations of Eqs. ( ). 

The integral along the imaginary axis contributes in the range 0 < Q 2 < oo. 

In numerical evaluations of the TPE correction, we use the analytical expressions for the 
master integrals of the poles contributions. The integration over the hyperangle 0 for the 
integral along the imaginary axis is performed numerically. 

The experimentally inaccessible region —Q 2 /4 < Q 2 < 0, corresponding with negative 
Q 2 values, is present in the Q 2 integration. For this relatively small region, we approximate 
the unpolarized proton structure functions by the following relations, F\ (w 2 ,— Q 2 ^j ~ 
Fi ^!E 2 ,g 2 j , F 2 (w 2 , — g 2 j ~ — F 2 (V 2 ,g 2 ), according to the approximation of Eqs. 
( ) and ( ). 
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FIG. 9: Q 2 dependence of the TPE correction 82 -, to e~p —> e~p for lab electron energy E = 
0.180 GeV (for which the kinematically allowed region is Q 2 < 0.094 GeV 2 ). The Feshbach term 
for pointlike particles, the elastic TPE contribution based on the box graph evaluation with dipole 
form factors [26], and the total TPE correction as the sum of elastic TPE and inelastic TPE 
are presented (upper panel). The inelastic contribution is compared with the leading logarithmic 
approximation of Eqs. ( ) and ( ) with A = 0.6 GeV and the result of Ref. [42] (lower panel). 

The experimental input for the proton structure functions is taken from the Christy-Bosted fit 
[48] . The vertical dashed lines restrict the region of validity of the expansion for the inelastic term 
Q 2 < (ME) /5 as follows from Fig. . 









36 




FIG. 10: Same as Fig. , but for the fixed value Q 2 = 0.05 GeV 2 as function of e in comparison 
with the empirical TPE fit using the data of Ref. [21] (A1 Collaboration, blue bands). 
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FIG. 11: Same as Fig. , but for the fixed value Q 2 = 0.25 GeV 2 . 


















